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ABSTRACT 



Context. The OH radical is a key species in the water chemistry network of star-forming regions, because its presence is tightly 
related to the formation and destruction of water Previous studies of the OH far-infrared emission from low- and intermediate-mass 
protostars suggest that the OH emission mainly originates from shocked gas and not from the quiescent protostellar envelopes. 
^ ' Aims. We aim to study the excitation of OH in embedded low- and intermediate-mass protostars, determine the influence of source 

, parameters on the strength of the emission, investigate the spatial extent of the OH emission, and further constrain its origin. 

. Methods. This paper presents OH observations from 23 low- and inteiTnediate-mass young stellar objects obtained with the PACS 

' integral field spectrometer on-board Herschel in the context of the "Water In Star-forming Regions with Herschel (WISH)" key 

program. Radiative transfer codes are used to model the OH excitation. 

Results. Most low-mass sources have compact OH emission (< 5000 AU scale), whereas the OH lines in most intermediate-mass 
sources are extended over the whole 47'.'0 x 47'.'0 PACS detector field-of-view (> 20000 AU). The strength of the OH emission is 
^ ' correlated with various source properties such as the bolometric luminosity and the envelope mass, but also with the [OI] and H2O 

emission. Rotational diagrams for sources with many OH lines show that the level populations of OH can be approximated by a 
CO I Boltzmann distribution with an excitation temperature at around 70 K. Radiative transfer models of spherically symmetric envelopes 

, cannot reproduce the OH emission fluxes nor their broad line widths, strongly suggesting an outflow origin. Slab excitation models 

'/^ ■ indicate that the observed excitation temperature can either be reached if the OH molecules are exposed to a strong far-infrared 

' continuum radiation field or if the gas temperature and density are sufficiently high. Using realistic source parameters and radiation 

fields, it is shown for the case of Ser SMMl that radiative pumping plays an important role in transitions arising from upper level 
, energies higher than 300 K. The compact emission in the low-mass sources and the required presence of a strong radiation field and/or 

a high density to excite the OH molecules points towards an origin in shocks in the inner envelope close to the protostar. 

Key words. Astrochemistiy — Stars: formation — ISM: molecules — ISM: jets and outflows 

X 

1. Introduction with varying temperature or radiation field, it is an excellent 

- ■ ■ probe of the physical conditio ns and dynamics in star-fo rming 

Oxygen is the most abundant element in the interstellar medium ^.^^^^^^ (Kristensen et al. 2012l :lvan Dishoeck et al.ll20Tlh . 
apart from hydrogen and helium. Many oxygen-bearing species, 

most importandy water and its precursors, have a veiy Um- A key connecting piece between atomic oxygen and water 

ited observability from the ground because of atmospheric con- is the hydroxyl radical (OH). In the high-temperature gas-phase 

straints. The Herschel Space Observatory dPilbratt et al .1120101) chemistry regime (T > 230 K), all available gas-phase oxygen is 

outperforms previous space-borne faciUties in sensitivity as well first driven into OH and then into water by the O H- H2 — » OH -1- H 

as spatial and spectral resolution. It is thus well suited to study and subsequent OH -H H2 — » H2O + H reactions. The impor- 

the water chemistry in young stellai" objects in more detail, tance of the backward reaction H2O -1- H — > OH + H2 de- 

Because water undergoes large gas-phase abundance variations pends on the atomic to molecular hydrogen ratio and there- 

fore on the local UV field or shock velocity. H2 is dissociated 

* Herschel is an ESA space observatory with science instruments in J- type shocks when the velocity exceeds ~ 25 km s"' (e.g. 

provided by European-led Principal Investigator consortia and with im- Ho llenbach & McKedl 19801) . OH is also a byproduct of the H2O 

portant participation from NASA. photo-dissociation process. Thus, OH is most abundant in re- 
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gions with physical conditions different from those favoring the 
formation and existence of water and it therefore provides a com- 
plementary view of oxygen in the gas phase. In addition to its 
importance in the oxygen and water chemistry, OH also con- 
tributes to the FIR cooling budget of warm gas in embedded 
YSOs (e.g., iN eufeld & Dalii^ [T989bl iKaufman & Neufeldl 
Il996t iNisini et al.l |200 2. Karska et al. subm.). The goal of 
the "Water In Star-Form ing regions with Herschel" (WISH, 
Ivan Dishoeck et ani201 Ih Herschel key program is to study the 
H2O and associated OH chemistry for a comprehensive picture 
of H2O during protostellar evolution. 

Detections of OH far-infrared (FIR) transitions from 
star-forming regions were made p reviously wi t h the 
Kuiper Airborne O bservatory (e.g. iMelnick et al.l Il987t 
Betz &Boreiko"l989'), t he Infrared Sp ace Obs ervatory (e.g. 
Ceccai-elli et al.. ,199a 



I, t ne inrrarea apace U Ds ervatory (e.g. 
I iGian nini et alj 20011; Larsson et al.l 



20021: iGoicoechea & Ce rnicharo ll2002l: l(3oicoechea et al.ll2004 



2006 ), Herschel (e.g. van Kempen et alJl2010al : I Wampfler et all 
20ld J2OII; Goicoechea et al. 20117 and recently SOFIA 



(ICsengerietal.ll2012t iWiesemeyer et al]l2012l) . Masers of OH 
hyperfine transitions are also commonly observed towards 
high-mass star-forming regions at cm wavelengths, but are not 
detected from low- and intermediate-mass protostars. Because 
this work focuses on low- and intermediate-mass young stellar 
objects (YSOs), we will not discus s OH rnaser emission (a 
detailed overview can be found in e.g. lElitzudll992h . 

From first Herschel results using t he Photodetecto r Array 
Camera and Sp ectrometer (PACS, iPoglitsch et"an l2010h . 
Ivan Kempen et alj j2010b,) found that the OH emission from the 
low-mass class I young stellar object (YSO) HH 46 is not spa- 
tially extended, in contrast to the H2O and high-7 CO emission 
from the same source. They speculated that at least parts of the 
OH emission could stem from a dissociative shock caused by the 
impact of the wind or jet on the dense inner envelope. Spectrally 
resolved Herschel observations of O H with the Heterodyne 
Instrument for the Far-Infrared (HIFI, Ide Graauw et al 1 120101) 
support an outflow scenario bas ed on the inferred broad lin e 
widths of more than 10 km s"^ d Wampfler et al.ll2010l 1201 ll) . 
Furthermore, analysis of OH PACS lines from a set of six low- 
and intermediate-mass protostars yielded similar e xcitation con- 
ditions of OH in all these sources (Wampfler et al. 2010). A ten- 
tative correlation of the OH line luminosities with the [OI] lumi- 
nosities as well as the bolometric luminosities of the protostars 
were found, indicating that the observed OH emission might 
originate from a dissociative shock. 

In this paper we present Herschel-PACS observations of OH 
in an extended sample of 23 low-and intermediate-mass YSOs. 
Our first goal is to test whether the tentative correlations from 
earlier work can be confirmed and to determine the influence of 
different source properties on the strength of the emission. The 
second goal is to study the OH excitation and the spatial extent 
of the OH emission in the target sources. A detailed analysis of 
the OH excitation and spatial extent is important to determine the 
origin of the OH emission and whether this is consistent with the 
picture from H2O observations. 

The paper is organized as follows: Section |2] describes the 
source sample, the observations, and the data reduction methods. 
In Sect.|3]we present the observational results. Section^Jcontains 
a discussion of spherically symmetric envelope models for OH 
(Sect. 14.11 ). the description and results of slab radiative transfer 
models to study the OH excitation in the outflow (Sect. 14.2b . as 
well as the discussion and data interpretation (Sect.|43). Finally, 
the conclusions are summarized in Sect. |5] 



2. Observations and data reduction 

Observations of 23 low- and intermediate-mass young stellar 
objects were carried out with the Photodetector Array Camera 
and Spectrometer (PACS, Poglitsch et al. 2010) on the Herschel 
Space Observatory. All observations were obtained within the 
"Water in Star-Forming Regions with Herschel" key program 
(WISH, van Dishoeck et al. 201 1) except for IRAS 12496-7650 
(DK Cha, ivan Kempen et alJ i2010al) . which was observed in 
the "Dust, Ice and Gas In Time (DIGIT)" key program (PI 
N. Evans). The coordinates and properties of the targets can be 
found in Table [T] and the data identity numbers (obsids), observ- 
ing modes, and the pipeline versions are given in Table IATI in 
the appendix. 

Two different observing modes are available for the PACS 
spectrometer, line and range spectroscopy. The "range spec- 
troscopy" mode provides a full scan of the 50 - 220 jum wave- 
length regime. The "line spectroscopy" mode covers only small 
windows around selected target lines, but generally at higher 
spectral sampling and sensitivity than range spectroscopy. The 
majority of the sources in our sample have been observed with 
the PACS line spectroscopy mode, targeting four main OH rota- 
tional doublets: ^ni/2(J = 1/2) ^ ^n3/2(7 = 3/2) at 79 yum, 
^n^/iiJ = 7/2 ^ 5/2) at 84 jum, ^UipM = 5/2 ^ 3/2) 
at 119 fim, and ^ni/2(J = 3/2 1/2) at 163 fim. The in- 
tegration time and the noise level for all sources observed in 
line spectroscopy mode is similar NGC 1333 IRAS 4A and 
NGC 1333 IRAS 4B have been observed in both fine and 
full range spectroscopy. Range spectroscopy only was used for 
NGC 1333 IRAS 2A, Ser SMM 1, and DK Cha. An iflustra- 
tion of the OH pure rotational transitions that are accessible with 
PACS on-board Herschel is provided in Fig.[T] An overview on 
the molecular data used in this work can be found in Table |2] 

The full spectral scan of DK Cha was presented previ- 
ously in 'van Kem pen et al.l (1201 Oal) and the line spectroscopy 
of HH 46 and NGC 7129 FIRS 2 in Fan Kempen et a l. { 201^ 
and Fich et al. (2010), respectively. These three spectra plus 
IRAS 15398, TMR 1, and NGC 1333 IR AS 2A were part of fli e 
sample analyzed in our previous work (IWampfler et al. I I20T0I) . 
We have now re-reduced all spectra with a newer version of 
the PACS pipeline and calibration files. The full spectral scans 
of NGC 1333 IRAS 4B and Ser SMMl are discussed ingreat 



detail inl Herczeg et al.l (1201 2h and IGoicoechea et aTl ( 1201 2|) re- 
spectively. 

The PACS spectrometer is an integral field spectrometer and 
operates simultaneously in a blue and a red channel. The detec- 
tor consist of 5 by 5 square spatial pixels ("spaxels") with a pixel 
size of 9'.'4. The Herschel half power beam width is smaller than 
a spaxel in the blue wavelength regime, but exceeds the spaxel 
size on the sky in the red part of the spectrum. The spatial reso- 
lution is therefore limited by the pixel size for short wavelengths 
and by the diffraction beam pattern at longer wavelengths. The 
spectral resolution depends on the grating order and varies from 
R = 3000-4000 at wavelengths ^ < 100//mto/? = 1000-2000 
for /I > 100 fim. All individual OH lines are therefore unresolved 
and at the lower resolutions, even blending of the OH doublets 
occurs, mostly at 79 fj.m and 119 jum. When doublet components 
are blended, they were assumed to be of equal strength in the 
analysis. 

The spectra were reduced with the H erschel interactive pro- 
cessing environment (HIPE, lOttI 1201 0^. versions 8 (for line 
scans) and 6 (for range scans). The wavelength grid was re- 
binned to four pixels per resolution element for line scans and 
two pixels per resolution element for range scans. The spectra 



2 



S. F. Wampfler et al.: OH survey with PACS 



Table 1. Coordinates, distance, bolometric luminosity, and envelope masses of the low-mass class 0, class I, and intermediate-mass 
protostars in the sample. 



Source 


RA 


Dec 


Class 


d 








[h m s] 


r ' "] 




[pc] 


[Lo] 


[Mo] 


NGC 1333 IRAS 2A 


03:28:55.6 


+31:14:37.1 





235" 


37.0" 


5.1" 


NGC 1333 IRAS 4 A 


03:29:10.5 


+31:13:30.9 





235" 


8.5* 


5.6" 


NGC 1333 IRAS 4B 


03:29:12.0 


+31:13:08.1 





235" 


4.9* 


3.0" 


L 1527 


04:39:53.9 


+26:03:09.8 





140'' 


1.8* 


0.9" 


CedllO IRS4 


11:06:47.0 


-77:22:32.4 





125'" 


1.0* 


0.2" 


IRAS 15398/B228 


15:43:01.3 


-34:09:15.0 





130" 


1.2* 


0.5" 


L 483 mm 


18:17:29.9 


-04:39:39.5 





200'' 


10.8* 


4.4" 


Ser SMMl 


18:29:49.8 


+01:15:20.5 





230^ 


31.7* 


16.1" 


Ser SMM3 


18:29:59.2 


+01:14:00.3 





230'^ 


11.0* 


3.2" 


L 723 mm 


19:17:53.7 


+ 19:12:20.0 





300'' 


3.9* 


1.3" 


L 1489 


04:04:43.0 


+26:18:57.0 


I 


140* 


3.7* 


0.2" 


TMR 1 


04:39:13.7 


+25:53:21.0 


I 


140* 


3.2* 


0.2" 


TMC lA 


04:39:34.9 


+25:41:45.0 


I 


140« 


2.7* 


0.2" 


TMC 1 


04:41:12.4 


+25:46:36.0 


I 


140* 


0.8* 


0.2" 


HH 46 


08:25:43.9 


-51:00:36.0 




450'' 


25.8* 


4.4" 


IRAS 12496/DK Cha 


12:53:17.2 


-77:07:10.6 




178' 


35.3* 


0.8" 


RN0 91 


16:34:29.3 


-15:47:01.4 




125^' 


2.4* 


0.5" 


AFGL 490 


03:27:38.4 


+58:47:08.0 


IM 


1000* 


2000 ' 


45 


NGC 2071 


05:47:04.4 


+00:21:49.0 


IM 


422" 


520 ° 


30 " 


Vela IRS 17 


08:46:34.7 


-43:54:30.5 


IM 


700' 


715 


6.4' 


Vela IRS 19 


08:48:48.5 


-45:32:29.0 


IM 


700' 


776 " 


3.5' 


NGC 7129 FIRS 2 


21:43:01.7 


+66:03:23.6 


IM 


1250' 


430 " 


50 " 


L1641 S3MMS1 


05:39:55.9 


-07:30:28.0 


IM 


465" 


70 


20.9^ 



References: " iHirota et aTI ( l2008h . * Values from lKristensen et al.l j2012h using PACS data from the WISH and DIGIT key programs (Karska et al. 
i n prep. Green e t al. in pr ep.), " Derived from DUSTY modeling of the sources as described in Kristensen et al. (2012), "' Andre et al. (2000), 

"'Knude & Hog' f 19981). ^ lEiroa et afl ( 120081) . Bas ed on VLB A measur e ments of a star thought to be associated with the S erpens cluster, 

bzib et al. (2010) inferred a distance of 415 pc. ^iK envonetai] J2008l), '' iHeathcote eTaP ( Il996h . ' IWhittet et aP (Il997 l). ^Ide Geus etJD (Il989l) . * 
Snell et al. ( 1984), ' Mozurkewich et al. ( 1986), "' Schreyer et al. (2002) d erived an envelop e mass of 40-50 Mp, " Wilson et al. ( 200 ^° 
Butner et aLU990.). Johnstone et al., C2001.). ''.Liseau et al., (il992). ' Giamiini et al.l (12003). 'iMassi et"^ ( ll999l) . ' iShevchenko & Yakubovl 
( I1989I) . " iJohnstone et al.l ( l2010h . ICrimieret al.l (120101) , "' IStanke et alj ( 120001) ' Ivan Kempen et alj ilOm 



Table 2. Molecular data from the LAMDA database (ISchoier et al ."2005') for the OH transitions detected with PACS. Wavelengths 
marked with a star were only observed in the full range scans. Frequencies are rest frequencies. A{B) = Ax 10*. 



Transition 




Wavelength 


Frequency 






Si. 


8t 


n, J, P 




[jjm] 


[GHz] 


[K] 


[s-'] 






3/2-3/2, 9/2-7/2, - - 


+ 


65.13* 


4602.9 


512.1 


1.276(+0) 


10 


8 


3/2-3/2, 9/2-7/2, + - 




65.28* 


4592.5 


510.9 


1.267(+0) 


10 


8 


1/2-1/2, 7/2-5/2, - - 


+ 


71.17* 


4212.3 


617.6 


1.014(+0) 


8 


6 


1/2-1/2, 7/2-5/2, + - 




71.22* 


4209.7 


617.9 


1.012(+0) 


8 


6 


1/2-3/2, 3/2-5/2, + - 




96.31* 


3114.0 


270.2 


9.270(-3) 


4 


6 


1/2-3/2, 3/2-5/2, - - 


+ 


96.37* 


3111.1 


269.8 


9.250(-3) 


4 


6 


1/2-1/2, 5/2-3/2, + - 




98.74* 


3036.3 


415.5 


3.530(-l) 


6 


4 


1/2-1/2, 5/2-3/2, - - 


+ 


98.76* 


3035.4 


415.9 


3.531(-1) 


6 


4 


1/2-3/2, 1/2-3/2, - - 


+ 


79.12 


3789.3 


181.9 


3.606(-2) 


2 


4 


1/2-3/2, 1/2-3/2, + - 




79.18 


3786.3 


181.7 


3.598(-2) 


2 


4 


3/2-3/2, 7/2-5/2, + - 




84.42 


3551.2 


291.2 


5.235(-l) 


8 


6 


3/2-3/2, 7/2-5/2, - - 


+ 


84.60 


3543.8 


290.5 


5.202(-l) 


8 


6 


3/2-3/2, 5/2-3/2, - - 


+ 


119.23 


2514.3 


120.7 


1.388(-1) 


6 


4 


3/2-3/2, 5/2-3/2, + - 




119.44 


2510.0 


120.5 


1.380(-1) 


6 


4 


1/2-1/2, 3/2-1/2, + - 




163.12 


1837.8 


270.1 


6.483(-2) 


4 


2 


1/2-1/2, 3/2-1/2, - - 


+ 


163.40 


1834.7 


269.8 


6.450(-2) 


4 


2 



were flat-fielded to improve the signal-to-noise ratio. The fluxes 
were normalized to the telescopic background and then cali- 
brated using measurements of Neptune as a reference. The rel- 
ative calibration uncertainty on the fluxes is currently estimated 
to be below 20%. The PACS spectrometer suff'ers from spectral 
leakage in the wavelength ranges 70 - 73 /vm, 98 - 105 yum. 



and 190 - 220 jum where the next higher grating order ranges 
52.5 - 54.5 fim, 65 - 70 yum, and 95 - 110 fj.m are superim- 
posed on the spectrum. The fluxes from lines in these wavelength 
bands might therefore be less reliable than in parts of the spectra 
that are not affected by spectral leakage. This applies to the OH 
71 jum and 98 jum doublets. 
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The lines were then subsequently analyzed in IDL using first 
or, if required, second order polynomials as baselines and the 
flux was measured by integrating over the line. The fraction of 
the point-spread function (PSF) seen by the central spaxel is 
wavelength-dependent, reaching from around 0.7 for a perfectly 
centered point source at 60-80 fim down to about 0.4 at 200 fim. 
A significant fraction of the flux might therefore fall onto neigh- 
boring spaxels and even more so if the source is off'-centered 
on the central spaxel. It is therefore important to investigate the 
flux distribution on the detector, which can be caused by spa- 
tially extended emission, the telescope PSF, or a combination 
of both. Extracting the flux from all 25 spaxels is not an opti- 
mal solution, because of contamination from nearby sources and 
because the signal-to-noise ratio drops if many spaxels without 
line emission but extra noise are added. This is particularly prob- 
lematic for weak lines or lines with little spatial extent. For a 
subset of our targets, where contamination from nearby sources 
occurs, we chose a set of spaxels that excludes the contribu- 
tion from close-by sources. This applies to NGC 1333 IRAS 4A, 
NGC 1333 IRAS 4B, Ser SMM3, and DK Cha. Spaxels excluded 
in the flux measurements of the WISH sources are marked in 
gray in the maps in the online appendixiBl For DK Cha, a 3 by 4 
set around the on-source spaxel was used. For all other sources, 
we used either the on-source spaxel, corrected for the spillover 
(L 1527, Ced 1 10 IRS 4, L 723, L 1489, RNO 91), 3 by 3 spaxels 
centered on the on-source spaxel (IRAS 15398, L 483, TMR 1, 
TMC lA, TMC 1, HH 46, and Ser SMMl), or flie full array (afl 
intermediate-mass sources). 
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Fig. 1. OH transitions accessible with Herschel/PACS. Wavelengths are given in units of microns. Transitions that were observed 
for the entire source sample are shown in red. Transitions targeted only in the fraction of sources observed in range scan mode are 
depicted in blue if detected from at least one source and in gray if undetected. 
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3. Results and Analysis 

3.1. Detected lines and spatial extent 

We detected at least one OH doublet in all 23 sources. The 
fluxes integrated over the emitting area, obtained as described in 
Sect. |2] can be found in Table [3] The doublet that is most often 
detected, in 21 out of 23 sources, is the ^113/2(7 = 7/2 — > 5/2) 
doublet at 84 yum, thanks to a combination of intrinsic strength 
and higher spectral resolution of the instrument at the shorter 
wavelengths. The component at 84.42 yum is however blended 
with CO(31-30) at 84.41 /im because the spectral resolution is 
around 0.037 fim. The sources in which the 84 jjm lines were 
not detected are NGC 1333 IRAS 2A and AFGL 490. Figures |2] 
and [3] present the OH spectra of the low-mass class and I 
sources. Fig. |4] the spectra of the intermediate-mass sources. In 
22 out of 23 sources at least one line was in emission, with the 
exception being NGC 1333 IRAS 2 A, where the only detec- 
tion i s the 119 fj.m doublet feature in absorption (IWampfler et al.l 
I2OIOI) . We mainly detected emission lines, but absorption is also 
found towards higher envelope masses. In the subsequent anal- 
ysis, we only considered the fluxes of lines that are purely in 
emission. Absorption in the 119 fim OH ^113/2 intra-ladder dou- 
blet transitions is observed from NGC 1333 IRAS 2A and in 
some spaxels of AFGL 490, Vela IRS 17, and Vela IRS 19, but 
there are also spaxels with emission from these targets. For even 
higher envelope masses like AFGL 490, the 79 fim OH cross- 
ladder transitions are also in absorption. This behavior can be 
explained by the fact that both lines are directly coupled to the 
ground rotational state of OH: because the first excited state is 
at Eup ~ 120 K, almost only the ground state is populated at 
temperatures below ~ 100 K. The 119 //m transitions couple the 
ground state to the first excited state and have a large Einstein A 
coefficient (~ 1.4 x 10"^ s"^). The 79 //m lines are cross-ladder 
transitions between the ground states of both rotational ladders. 
Cross-ladder transitions generally have much smaller Einstein A 
coefficients (~ 3.6 x 10"^ s"' for the 79 fim lines) than the intra- 
ladder transitions. Therefore, absorption at 79 fim does not occur 
as easily as for the 119 fim doublet. 
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Table 3. OH fluxes used throughout the paper in units of 10 W m ^ with Icr errors from the flux extraction. Calibration errors are not included. Values in round brackets are 
detected below the 5cr but above the So" level, those in square brackets are above 2cr but below 3cr. 



Source 



OH transition [/um] 





79.12 


79.18 


84.42 


84.60 


119.23 


1 19.44 


163.12 


163.40 


NGC 1333 IRAS 2A (scan) 


- 


- 


- 


- 


abs. 


abs. 


- 


- 


NGC 1333 IRAS 4A 


(0.98 + 0.23) 


0.99 ± 0.23 


2.40 ± 0.16 


0.84 ±0.14 


[0.41 ± 0.15] 


1.51 ± 0.15 


0.64 + 0.07 


0.53 + 0.07 


NGC 1333 IRAS 4B 


2.08 ±0.16 


1.30 ±0.14 


4.88 ±0.19 


2.50 ±0.17 


2.23 ± 0.09 


2.77 ± 0.08 


0.82 + 0.07 


(0.58 + 0.08) 


L 1527 


0.33 ± 0.05 


0.27 ± 0.05 


0.50 ± 0.05 


0.49 ± 0.06 


0.60 ± 0.03 


0.71 ± 0.03 


[0.07 + 0.03] 


(0.12 + 0.03) 


CedllO IRS4 


0.31 ± 0.06 


(0.23 ± 0.06) 


0.53 ± 0.06 


0.47 ± 0.06 


(0.16 ±0.04) 


0.25 ± 0.06 


0.16 + 0.03 


- 


IRAS 15398/B228 


0.63 ± 0.14 


0.65 ± 0.14 


1.18 ± 0.13 


0.58 ±0.12 


0.73 ± 0.09 


0.69 ± 0.08 


- 


- 


L 483 mm 


1.17 ± 0.18 


1.15 ± 0.18 


1.46 ± 0.14 


0.88 ± 0.14 


0.36 ± 0.06 


0.65 ± 0.07 


0.31+0.04 


(0.17 + 0.04) 


SerSMMl (scan) 


9.97 ± 0.82 


blend 


12.5 ± 0.43 


9.29 ± 0.43 


8.75 ± 0.71 


8.18 + 0.71 


2.76 + 0.30 


1.92 + 0.33 


Ser SMM3 


(1.46 ± 0.32) 


1.57 ± 0.31 


2.01 ± 0.30 


1.38 ± 0.29 


(abs.) 








L723 mm 


- 


[0.17 + 0.06] 


0.61 ± 0.06 


0.41 + 0.05 


(0.12 + 0.03) 


(0.13 + 0.04) 


- 


(0.14 + 0.03) 


L 1489 


0.79 ± 0.07 


0.77 ± 0.06 


1.03 ±0.07 


0.94 ± 0.08 


0.52 ± 0.03 


0.67 + 0.02 


(0.14 + 0.04) 


(0.15 + 0.04) 


TMR 1 


(0.72 ±0.16) 


1.17±0.16 


1.35 ±0.16 


1.12 ±0.16 


(0.35 ± 0.08) 


0.47 ± 0.08 


0.30 + 0.05 


0.28 + 0.05 


TMC lA 


- 


- 


0.77 ± 0.13 


0.66 ± 0.13 


(abs.) 


(abs.) 


- 


- 


TMC 1 


1.00 ±0.08 


0.89 ± 0.08 


1.40 ± 0.14 


1.11 ± 0.14 


0.45 ± 0.07 


0.67 ± 0.06 


(0.18 + 0.06) 


- 


HH 46 


0.69 ±0.17 


0.73 ± 0.18 


0.84 ±0.17 


1.00 ± 0.18 


0.37 + 0.06 


0.58 + 0.06 






DK CHA 


(2.29 ± 0.68) 


blend 


4.82 ± 0.44 


2.88 ± 0.40 










RN0 91 


0.61+0.11 


blend 


(0.36 ± 0.07) 


(0.40 ± 0.08) 




[0.08 + 0.03] 






AFGL490 


mix. 


mix. 






abs. 


abs. 


1.19 + 0.08 


0.90 + 0.08 


NGC 2071 


30.96 + 0.73 


30.47 + 0.79 


42.54 + 0.81 


33.68 ± 0.75 


mix. 


mix. 


12.62 + 0.39 


8.33 + 0.37 


Vela IRS 17 


2.60 ± 0.23 


2.50 ± 0.24 


3.12 ±0.27 


2.92 ± 0.24 


mix. 


mix. 


0.95 + 0.13 


(0.51 + 0.12) 


Vela IRS 19 


3.08 ± 0.44 


2.62 ± 0.37 


4.99 ± 0.49 


2.20 ± 0.44 


mix. 


mix. 


1.27 + 0.07 


1.01+0.08 


NGC 7129 FIRS 2 


2.60 ±0.33 


(1.41 ± 0.36) 


1.86 ±0.30 


(1.30 ±0.26) 


1.05 + 0.11 


1.20 + 0.13 


0.63 + 0.08 


(0.42 + 0.08) 


L1641 S3MMS1 


2.49 ± 0.43 


2.93 ± 0.45 


4.28 ±0.31 


2.87 ± 0.30 


mix. 


mix. 


0.87 + 0.08 


0.93 + 0.08 
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J I Ll_l Ll_l I I I I LJ I L 



79.1 79.2 84.4 84.6 119.2 119.6 163.2 163.6 

Wavelength A (/^m) 

Fig. 2. PACS line scan OH spectra of the low-mass class young stellar objects in our sample. The x-axis is wavelengths in yum, 
the y-axis is continuum subtracted flux density in Jy (plus a constant off'set). The red dashed lines indicate the rest frequencies of 
the OH transitions. The blue dotted line represents the rest frequency of CO(31-30), which is blended with OH at 84.42 fim. The 
sampling of range scans (Ser SMMl and NGC 1333 IRAS 2A) is different from line scans. 
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Fig. 3. The same as Fig.|2]but for the low-mass class I YSOs. 
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Fig. 4. The same as Fig.|2]but for the intermediate-mass protostai^s. 
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Fig. 5. Map of the OH 84 //m (red) and [OI] 63 fj.m emission 
(black) from NGC 2071, illustrating the similar spatial extent of 
the two transitions. The x-axis is velocity in km s"', the y-axis 
continuum subtracted and normalized flux density. All spaxels 
were normalized with respect to the spaxel that contains the peak 
of the continuum emission at 63 jum. Note that the peak of the 
continuum emission at 84 fim falls onto the central spaxel, al- 
though the two observations were carried out consecutively. 



Several sources show OH emission that is extended beyond 
what would be expected from a point source folded with the 
point-spread function of the telescope or leakage onto neigh- 
boring spaxels. The line emission is usually extended along the 
outflow direction (see Karska et al. subm.) and strongly corre- 
lated with the spatial extent of the atomic oxygen transitions, as 
illustrated in Fig.|5]and further discussed in Karska et al. (subm). 

Maps of the OH 79,84,119, and 163 yum transitions for 
all WISH sources can be found in the appendix. An example 
of compact OH emission at 84 /im from the low-mass class I 
YSO L 1489 is shown in Fig.|6] Despite their location at further 
distances than the low-mass objects, all six intermediate-mass 
sources show signatures of extended OH emission. The largest 
spatial extent is seen in the 119 yum ground state lines, which 
can be excited most easily. Figure |7]presents the spaxel map of 
the 119 yum line from AFGL 490. The absorption is strongest to- 
wards the YSO and extended over an area of ~ 25" (25000 AU) 
around the central position in the direction perpendicular to 
the outflow. The spatial extent of the emission is comparable 
to the 20000 AU by 6000 AU envelope structure discussed in 
ISchrever et al.l (l2006h . The OH absorption changes into weak 
emission along the outflow direction. For the full sample we 
provide an overview whether the flux observed outside the on- 
source spaxel is consistent or inconsistent with the spillover fac- 
tor in Table g] 



3.2. OH emission line flux ratios 

Comparison of the OH 84 jjm line luminosities (flux corrected 
for the source distance) between the different source types shows 
that the intermediate-mass sources in our sample have higher OH 
luminosities than the low-mass sources by about two orders of 
magnitude on average. Among the low-mass sources, the class 
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Fig. 6. Map of the compact OH 84 jim emission from L 1489. 
The X-axis is wavelength in yum, the y-axis continuum subtracted 
and normalized flux density. All spaxels were normalized with 
respect to the central one. The red dashed lines indicate the 
rest frequencies of the OH transitions, the blue dashed line the 
CO(31-30). 
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Fig. 7. Map of spatially extended OH 119 fim absorption and 
emission from AFGL 490. The x-axis is wavelength in fim, the 
y-axis continuum subtracted and scaled flux density. The labels 
indicate the multiplicative scaling factors for each spaxel unless 
they are unity. 



sources are on average about a factor of two more luminous in 

OH than the class I s ources. 

In earlier work (IWampfler et al 1 1201 Oh . we found that the 
line flux ratios among the sources in the sample were relatively 
constant. We have therefore also calculated the 79 yum/84 fim, 
79 fim/ll9 fim, 79 fim/163 fim, 84 yum/1 19 fim, 84 fim/163 fim, 
and 119 fim/163 fim line flux ratios for the extended sample, i.e. 
the sources in Tab. [1] from which emission was detected. The 
fluxes of doublets are added except for the 84 yum doublet, be- 
cause the 84.42 yum is blended with CO(31-30). Instead we use 
twice the flux of the 84.60 yum component assuming that both 
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Table 4. Description whether the flux measured in the on-source spaxel and the surrounding three by three spaxels are consistent 
with the pure spillover factor (0.70 at 79 jum, 0.69 at 84 jum, 0.62 at 1 19 fim, and 0.49 at 163 fim). If the fraction of flux measured 
in the on-source spaxel compared to the three by three central spaxels is smaller than the spillover value, it is an indication that the 
OH line is spatially extended, as summarized in the last column. 



Source 


79 yum 


84 yum 


119yum 


163 /jm 


spatially extended 


NGC 1333 IRAS 2 A 






no (abs.) 






NGC 1333 IRAS 4 A 




no 


no 




yes 


NGC 1333 IRAS 4B 


no 


no 


no 


no 


yes 


L 1527 


no 


VPS 

yes 


yes 


yes 


Cno'l 


CedllO IRS4 


yes 


( ves^ 


yes 


yes 


no 


IRAS 15398/B228 


no 


no 


no 


no 


yes 


L 483 mm 


(no) 


ves 


yes 




(no) 


SerSMMl 


no 


no 


no 


no 


yes 


Ser SMM3 


no 


no 






yes 


L 723 mm 




yes 


yes 




no 


L 1489 


yes 


VPS 
yes 


yes 


yes 


no 


TMR 1 


no 


no 


no 


yes 


(mispointed) 


TMC lA 


no 


yes 




yes 


(no) 


TMC 1 


no 


no 




yes 


(yes) 


HH 46 


no 


yes 


no 




no 


DK Cha 




no 






(yes) 


RN0 91 


yes 








no 


AFGL 490 


no (abs.) 


no (abs.) 


no (abs.) 


no 


yes 


NGC 2071 


no 


no 


no (abs./em.) 


no 


yes 


NGC 7129 FIRS 2 


no 


yes 


no 




yes 


Vela IRS 17 


no 


no 


no (abs./em.) 


no 


yes 


Vela IRS 19 


no 


no 


no (abs./em.) 


no 


yes 


L1641 S3MMS1 


no 


yes 


no 


yes 


yes 



tected in the full spectral scan of NGC 1333 IRAS 2 A and only 
very few line detectio ns are available for NGC 1333 IRAS 4 A 
and DK Cha (see also Ivan Kempen et alj|20103) . so that the ro- 
tational diagram for these sources is very sparsely populated 
and a fit of the rotational temperature is therefore not feasible. 
For Ser SMMl, the corresponding OH column density would 
be A^oH = 1.0 X lO''* cm"^ assuming a source size of 20" and 
using an inte rpolated value for th e partition function Q from the 
JPL catalogFl dPickett et al.ll998l) . However, the ^03/2 points fall 
below the fit, suggesting that these transitions are either opti- 
cally thick, which is well possible because the ^113/2 ladder con- 
tains the ground rotational state, or that the two rotational ladders 
might have a different rotational temperature. 

3.4. Dependence of OH luminosity on source properties 

Understanding which source parameters influence or even deter- 
mine the strength of the OH emission and the excitation of the 
different transitions is important to constrain the origin of the 
emission. We therefore test whether the OH line luminosity cor- 
relates with various envelope parameters and the emission from 
other molecular and atomic species. Because the different OH 
transitions are fairly well correlated (Sect. 13.2) . we restrict the 
analysis to the 84 fim luminosity, where we have most detec- 
tions. 

In Fig.|9]and Fig.[TO]the dependence of the OH line luminos- 
ity on the bolometric source luminosity and the envelope mass 
are shown, respecti vely. We use new values for Lboi that were de- 
rived by Kristense n et al.l (1201 2h based on additional continuum 
values from PACS observations, which are presented in Karska 
et al. (subm.). The envelope masses for the low-mass sources 
are calculated from spherical models based on continuum radia- 
tive transfer and include all material at temperatures above 10 K. 

' http://spec.jpl.nasa.gov 



Table 5. Observed line flux ratios (see text for details). 



Ratio 


Minimum 


Maximum 


Mean 


Median 


79/ 84 


0.40 


1.54 


0.90 


0.85 


79/119 


0.46 


2.30 


1.32 


1.31 


79/163 


1.68 


5.38 


3.22 


3.01 


84/119 


0.75 


3.28 


1.65 


1.58 


84/163 


1.44 


6.48 


3.44 


3.57 


119/163 


1.41 


4.10 


2.66 


2.14 



components are of similar strength. Cases where only one dou- 
blet component was clearly detected were not considered. The 
results are listed in Table |5] 

Again we find that the line ratios remain relatively constant 
within less than a factor of four around their mean values over 
the whole luminosity, mass, and age range of several orders of 
magnitude spanned by the sample. The excitation of OH is there- 
fore likely to be similar in all sources, indicating that either the 
OH emission stems from gas at similar physical conditions or 
that the ratios remain stable over a spread of parameter values. 
The latter possibility is supported by the models (cf. Sect. |4]l, but 
does not exclude the first option. 

3.3. OH rotational temperature 

The full spectral scans cover the OH transitions from about 55 - 
200 jum (Eap ~ 120 - 875 K) and therefore allow us to study the 
excitation conditions by means of rotational diagrams. Figure [8] 
presents the diagram for Ser SMMl . The derived rotational tem- 
perature is Trot ~ 72 ± 8 K if all the transitions shown on the 
plot except the 119 yum (optically thick), 84.42 yum (blended with 
CO), and 98 |/m (in leaking region) are included. The rotational 
diagra m for NGC 1333 IRAS 4B can be found in lHerczeg et al.l 
(120 12h . giving T^rot ^ 60 ±15 IC. No emission lines were de- 
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Fig. 8. Rotational diagram for the OH lines measured from 
Ser SMMl. Transitions belonging to the '^Tli/2 ladder are shown 
in red, those from the ^03/2 ladder in blue, and cross-ladder tran- 
sitions in green. The error bars indicate the 3 cr errors on the flux 
measurement (calibration errors are not included). 
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Fig. 9. OH luminosity of the 84 fim transition vs. bolometric lu- 
minosity of the source. Low-mass class sources are shown as 
red diamonds, class I sources as blue triangles, and intermediate- 
mass protostars as green circles with 3 cr error bars from the flux 
determination (does not include the calibration uncertainty). The 
dashed line is a linear fit to the data. 



For the intermediate-mass sources, we use literature values and 
the method that was used to determine the envelope mass may 
therefore be different. 

The Pearson correlation coefficient, defined as 
Px,Y = cov{X,Y)/[cr{X) ■ o-(Y)], i.e. the covariance of two 
random variables X and Y divided by their standard deviations, 
for LoH with Lboi is 0.93, including all 21 sources where the 
OH 84 fj.m line was detected. An overview on all obtained 
correlation coefficients and their corresponding significance 
levels is given in Table|6] Values of p close to 1 and -1 describe 
a tight correlation or anticorrelation of the parameters X and 
Y, respectively, while values close to indicate that X and Y 
are uncorrected. At what level of p a correlation is considered 
to be significant, depends on the number of data points. The 
significance of p can be expressed as a probability value, the 
significance level p, and tells how likely the actual observed or a 
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Fig. 10. OH 84 fim luminosity vs. envelope mass of the source. 
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Fig. 11. Bolometric luminosity vs. envelope mass of the sources. 

Table 6. Pearson's correlation coefficients p, number of sources 
A^, corresponding significance levels p, and number of standard 
deviations cr for the OH 84.6 yum line luminosity or flux and 
various envelope parameters. 



Parameter 


P 


A' 


P 


cr 




0.93 


21 


< 10 ' 


4.17 


-A^eiiv 


0.84 


21 


< 10-3 


3.77 


Tbo\ 


-0.23 


16 


0.20 


0.87 




0.46 


15 


0.04 


1.71 




-0.21 


21 


0.18 


0.95 


100 K model radius 


0.86 


16 


< 10-3 


3.33 


"1000 AU 


0.47 


16 


0.03 


1.82 


Ht column density 


0.30 


16 


0.26 


1.16 


[01] 63 fim 


0.71 


21 


< 10-3 


3.16 


[01] 145 urn 


0.79 


18 


< 10-3 


3.27 


H2O 89.99 i^m 


0.83 


17 


< 10-3 


3.33 



more extreme value of p is found under the assumption that the 
null hypothesis is true, i.e. that there is no relationship between 
the parameters. Another option is to express the significance 
in terms of number of standard deviations cr, calculated from 
p yN-l, where is the number of sources. 

The correlation coefficient for Lqh with Menv is smaller 
(p - 0.84). As illustrated by Fig. (TT] the bolometric lumi- 
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Fig. 12. OH 84 fim luminosity vs. bolometric temperature of the Fig. 13. OH 84 /zm luminosity vs. outflow force, 
sources. 



nosity and the envelope mass are not independent properties 
of the sourc es. This is c onsistent with iBontemps et al.l (1 19961 
Fig. 1) and Andre et al.l (l200ft Fig. 6b), who found that Lboi 
and Man are strongly correlated for the class sources and 
that class 1 sources move down in the diagram with time along 
the evolutionary tracks. From an Lqh - A^env correlation one 
would also expect L qh to be correlated with Tboi (Fig. 8 of 
lJ0rgensen et al.ll2b02h . Such a correlation is however not found 
for our low-mass source sample, as illustrated by Fig. [12] in- 
dicating that the underlying Lboi - Me,„ correlation might cre- 
ate the observed Lqh - A^env trend. A correlation between two 
variables can be caused by an underlying correlation of the two 
variables with a third one. A method to measure such effects 
is the concept of the partial correlation coefficient, defined as 

Pixjyz = (Pxj - Px,z ■ Py,z)I ^ (1 - pj^z) ' ~ Py.z) (^^^ ^1^° 
appendix A from lMarseille et ani2010l) . The resulting value for 
the partial correlation coefficient of Lqh with Menv under the in- 
fluence of Lboi is 0.56, i.e. significantly reduced, thus indicating 
that the trend might indeed be based on an underlying relation 
between mass and luminosity. Class 1 YSOs that lie above the 
LoH - ^env Correlation in Fig. [10] are also the ones that do not 
follow the Lboi-A^env relation (Fig. [TOl l. in particular DK Cha. 
These sources have a higher Lqh than what would be expected 
from their mass and a lower value than what would be expected 
from their bolometric luminosity. They are found in the regio n 
of more evolved class 1 sources (see Fig. 6 in lAndre et al.ll2000l) . 
i.e. they have a significantly lower envelope mass at a given bolo- 
metric luminosity than less evolved sources like e.g. HH 46, and 
are in the transitional stage to class 11 where disk emission may 
contribute as well. 

We find that Lqh (the luminosity of the OH 84 jum tran- 
sitions) seems to be well correlated with luminosity for both 
evolutionary stages and that class and 1 s o urces lie on 
the same straight line. In contrast, iNisini et al ] (120021) found 
from ISO data that class I sources fall systematically below 
class sources in their plot of the total far-infrared luminosity 
LpiR = ^^[01] + Leo + i'HoO + Lqh versus Lboi- Furthermore, they 
concluded that class and 1 sources fall on to the same straigh t 
line in the Lfir - Me„v plot, but in Fig. 6 of iNisini et al.l (l2002h . 
class 1 sources with Me„v < 10"' tend to branch off as well. 

A correl ation of Lqh with the o utflow momentum flux ("out- 
flow force" iBontemps et al.lll996b Fqo is not significant from 
our Fig.[T3] The values for Fco are taken from the literature and 




J 10-^ 



Fig. 14. OH 84 /zm luminosity vs. luminosity to the 0.6 divided 
by the envelope mass of the sources (evolutionary tracer). 



were not derived in a consistent way. New estimates of Fqo are 
currently work in progress (Yildiz et al. in prep.). 

Furthermore, Lqh seems to decrease with evolution for 
class sources, but not for class 1, as can be seen from Fig. [14] 
showing Lqh plotted against l}^^^y| Mf.,nw, a quantity t hat w as pro- 
posed as an evolutionary tracer bv lBontemps et al.l(ll996l) . 

The different behavior of the class and I sources could be 
caused by their different geometry. If the OH was radiatively 
pumped, then the location of the warm dust would be an im- 
portant parameter Warm dust at 100 K (cf. Sect. 13. 3t is ex- 
pected from two physical components in young stellar objects, 
the inner envelope and the outflow walls. Class and class I 
sources differ in the geometrical alignment of their warm dust 
components: In class sources, high-energetic radiation from 
the protostar and accretion processes is reprocessed to a far- 
infrared field by the dust in the inner envelope, which repre- 
sents a small solid angle. Thus the OH luminosity for sources 
where this component dominates should be mass dependent. In 
class 1 sources, the warm dust from the outflow walls can di- 
rectly irradiate the local OH molecules and therefore the OH 
luminosity for these sources should depend on the bolometric 
luminosity. We extracted the radius at which the temperature 
reaches 100 K from the spherically symmetric envelope mod- 
els to test this hypothesis. Figure [TS] shows the OH luminosity 



14 



S. F. Wampfler et al.: OH survey with PACS 



s 

x 



-3 



10"^ - 



p - 0.86 



♦ 



10 



100 



100 K Radius (AU) 



-14 



10 



7 io-'« 



S 10 



-16 



10 



-17 



p = 0.71 



o 



63 /xm 



■ I ■ 



p = 0.79 o 



145 fxrci 



LUI ■ I 



10 



10 



-14 



Fig. 15. Correlation of OH 84 ^im luminosity with the 100 K Fig. 17. Correlation of OH 84 fxm and [OI] 63 /im and 145 jum 
radius of the spherical source models. fluxes. 




10 



-14 



1000 AU Density (cm-^) 



o 



X 10-16 



10 



-17 




F„2o (W 



Fig. 16. Correlation of OH 84 luminosity with the H2 density Fig. 18. Correlation of OH 84 //m and H2O 89.99 //m fluxes, 
at 1000 AU from the spherical source models. 



plotted against th e 100 K radius from th e spherically symmetric 
envelope models ( Kristensen et alj2012l) and there seems to be a 
trend of increasing OH luminosity when the size of the inner hot 
envelope T > 100 K is larger The OH luminosity does not seem 
to depend significantly on the density at 1000 AU as shown in 
Fig- EH A correlation with density would not necessarily imply 
that collisions are the dominant excitation mechanism anyway 
as the dust mass also increases with density. 

3.5. Correlation of the OH flux with other species 

OH is tightly related to the formation and destruction processes 
of H2O and a connecting piece between H2O and atomic oxy- 
gen in the chemistry. We therefore test whether the OH emis- 
sion is correlated with the [OI] and H2O emission. The [OI] and 
H2O fluxes for the low-mass protostars are tabulated in Karska 
et al. (subm.) and those of the intermediate-mass protostars in 
Table lA.2l in the appendix. Fig.[T7]shows the OH 84 fim flux plot- 
ted versus the [OI] ^Pi - ^P2 and ^Po - ^^Pi fluxes at 63 fim (left 
panel) and 145 jum (right panel), respectively. The Pearson corre- 
lation coefficients for OH 84 fj.m vs. [OI] 63 yum and 145 jum are 
0.71 and 0.79, respectively, corresponding to a three sigma re- 
sult. The [OI] 63 fim could be optically thick or more affected by 



contamination than the 145 /im, which would explain the slightly 
lower correlation coeflicient for OH 84 fj.m vs. [OI] 63 fim com- 
pared to [OI] 145 fim. The OH and [OI] emission thus seem to 
be correlated both spatially and in total flux, hinting at a com- 
mon origin, most likely from gas associated with the outflow, 
because in the cases of non-compact emission it is usually found 
to be extended along the outflow direction. 

The correlation of the OH 84 fim flux with P-H2O 32 2 - 2i 1 
(A ^ 89.99 yum, E^p = 296.8 K) is shown in Fig. [ft] The 
P-H2O 32,2 - 2ij was chosen because it has a very similar up- 
per level energy to the OH 84 fim transition (^up - 290.5 K) 
and a similar wavelength, which strongly reduces the influence 
of instrumental effects, as the PSF is very similar. The fluxes 
were measured on the same spaxel sets. The correlation of the 
OH 84 jum flux with P-H2O 32,2 - 2 1,1 is significant at more 
than 3cr with p = 0.83. The same plot for 0-H2O 2i,2 - lo,i 
{A X 179.53 //m, ^up = 114.4 K) can be found in Karska et 
al. (subm.). The 0-H2O 2i,2 - lo,i line has a lower upper level 
energy and was observed in a larger beam, but was detected in 
more sources. Again, there is a statistical correlation between the 
fluxes of the OH and H2O lines at similar energy, but the H2O 
and OH spatial extents can differ (Karska et al. subm.). 

The ratio of the two fluxes versus the bolometric luminosity 
is shown in Fig. [19] The low-mass class I sources lie in the up- 
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per part of the plot as they have OH 84 yum to H2O 89 yum ratios 
above 2-3, whereas the class sources have a slightly lower 
ratio on average. On the other hand, there is no significant trend 
with bolometric luminosity (p = -0.34). The ratio seems to de- 
crease with envelope mass, as illustrated by Fig. |20] and to in- 
crease with bolometric temperature (Fig. |211 . but there is a large 
spread. Fig. |22] shows that an increase in the OH 84 yum to H2O 
89 fim ratio could be an evolutionary effect. The increase in OH 
compared to H2O could for instance be caused by the envelope 
becoming more tenuous in the more evolved stages, allowing 
high-energy photons from the protostar to penetrate further and 
dissociate the H2O. The intermediate-mass sources are missing 
in Fig.|2Tl because a consistent set of bolometric temperatures is 

not available. 

lOoicoechea et al.l (1201 ll) find that the OH 84 yum emission 
is not spatially correlated with 0-H2O 3o,3 - 2] 2 in the Orion 
Bar PDR, but with high-7 CO and CH+. We cannot test such a 
spatial correlation here because the fluxes drop rapidly outside 
the on-source spaxels in our low-mass source sample. 



4. Excitation and origin of OH 

Where does the OH emission arise? The environment of low- 
mass protostars is known to consist of dififerent physical com- 
ponents, including the collapsing envelope, shocks, outflow cav- 
ities heated by UV radiation and entrained outflo w gas, all of 
which are contained in the Herschel beams (see IVisser et al.l 
2012, for an overview). To derive physical properties like the 
density and temperature of the emitting gas, the observed fluxes 
and line ratios are compared to the results from radiative transfer 
models. The obtained physical conditions can help in constrain- 
ing the origin of the emission in combination with information 
about the spatial distribution of the emission. Furthermore, ra- 
diative transfer models permit to study what mechanism domi- 
nates the excitation under given physical conditions. 

The lowest rotational transitions of OH lie in the far-infrared 
and therefore significantly interact with the dust continuum field, 
which peaks in the same wavelength regime for dust temper- 
atures typical of embedded YSOs. Moreover, the properties of 
OH such as the large dipole moment and large rotational con- 
stant make radiative pumping an important excitation process. 
Therefore, the excitation of OH should be modeled using a code 
that includes radiative effects. 
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Before describing the models, we recall the clues on the 
origin of the OH emission that come from spectrally resolved 
OH line profiles. Wampfler et al. (2010) present upper limits 
on the OH 163 //m (1837 GHz) line intensities using HlFl for 
two low-mass sources which, when combined with measured 
PACS fluxes for the same lines, imply FWHM line widths of 
at least 10 km s"' . For a third low-mass source, Ser SMMl, a re- 
cent unpublished HlFl spectrum shows a detection of the broad 
(FWHM ^ 20 km s ') component, but no sign o f a narrow com- 
ponen t (Wampfler et al., in prep.). In contrast, IWampfler et al.l 
Hon]) found a narrow component (FWHM =s 4-5 km s ^) in the 
OH 1837 GHz spectrum from the high-mass protostar W3 IRS 5 
which can be attributed to the quiescent envelope, on top of a 
broad outflow component with FWHM * 20 km s ' . Therefore, 
the OH line profiles from low-mass YSOs seem to be dominated 
by the outflow contribution, whereas those of high-mass sources 
can contain both an outflow and envelope component. 



4.1. Spherical envelope models 

To further quantify any contribution from the protostel- 
lar envelop e, spherically symmetri c source models such as 
derived by fcistensen et al.l (1201 2|) are commonly used in 
combination wit h a non-LTE line radiat i ve tra nsfer code 
like RATRAN (Hogerheiide & van der Tag |2000) or LIME 
(Brinch & Hogerheiide 2010) to model the line emission. 

We ran RATRAN calculations in spherical symmetry for a 
few low-mass sources, including Ser SMMl, using the same 
RATRAN code as for W3 IRS 5 (see Sect. 14X11 for molecu- 
lar data input). The OH fractional abundances considered range 
from xoH = 10"'' to 10"^. Note that the source structure on small 
scales is often not well constrained by the single-dish contin- 
uum observations and that the s pherical symmetry star ts to break 
down inside about 500 AU (e.g.'j0rge nsen et al.l2005 !). Thus this 
type of model is not ideally suited for lines with a high criti- 
cal density, where the bulk of the emission usually arises from 
the innermost parts of low-mass protostellar envelopes. The in- 
ner envelope is also the area where the excited states of OH are 
mainly populated, and thus the model results should only be re- 
garded as a rough order of magnitude estimate. 

It is clear that these spherically symmetric source mod- 
els cannot reproduce the observed OH emission in low-mass 
sources for a variety of reasons. First, the synthetic spectra from 
the model fail to reproduce the observed broad lines (FWHM 
of ~ 20 km s '). The model lines are typically much narrower 
(FWHM of a few km s ') for any reasonable Doppler-fe pa- 
rameter that fits emission from other molecules. Therefore, the 
line width measured from the HIFI spectra is inconsistent with 
an envelope-only or envelope-dominated scenario for the low- 
mass sources observed with HIFI. Moreover, the 79, 84, and 
119 jum lines from the RATRAN models are always in absorp- 
tion, while the detected PACS lines from the low-mass sources 
are in emission, providing another indication that the bulk of the 
emission does not arise from the envelope. Although a minor 
envelope contribution to the 163 fim lines cannot be excluded, 
the bulk of the emission is likely associated with the outflow. 
For the intermediate-mass sources, absorption is observed at 
119 yum in at least one spaxel (except for NGC 7129 FIRS 2), 
suggesting that a potential envelope contribution to the spectra 
of intermediate-mass sources cannot be excluded, in particular 
in the lowest rotational transitions. 



4.2. Outflow model 

4.2.1. Slab model Description 

To model OH emission from outflows or shocks, a simple 
slab geometry is appropriate. We use a radiative transfer code 
based o n the escape probability formalism as described in 
Takahas hi et al. l (fT98l . which includes radiative pumping. The 
physical conditions such as density and temperature and the ex- 
citation of the molecule are assumed to be constant through- 
out the region. The free parameters of the model are the ki- 
netic temperature of the gas Tgas, the temperature of the dust 
Tdust, the density of the collision partner «, the molecular col- 
umn density of OH A^oh, and the dust column density A^dust- 
The line profile function has a rectangular shape with a width 
of Ay = yfTT/il Vln 2)Av where Av is the full width at half max- 
imum (FWHM) of a Gaussian line. Here, Av was assumed to 
be 10 km s"' for all models. As discussed above, this FWHM 
is motivate d by the non-det ection of OH from low-mass sources 
with HIFI llWamnfleretanilomi that lead to the conclusion that 
the line width must be broader than 10 km s ' . 

Our slab co de is similar to the "RADEX" code 

(Ivan derTaketal.1 l2007h . but includes an active dust con- 
tinuum, which is capable of both absorbing and emitting 
radiation. Results from our code in the limiting case where no 
dust is present agree very well with the RADEX resu l ts. An 
earlier OH excitation study by Off er & van Dishoeckl (Il992h 
included a dust continuum radiation field for the radiative 
excitation but not the absorption of line photons by dust 
grains. Furthermore, different dust properties were used in their 
calculations. Therefore, the comparability of the two models is 
limited. 

We use the molecular data file from the Leiden atomic and 
molecular database (LAMDA, Schoi er et al. 2005) for OH with- 
out hyperfine structure, because the resolution of PACS does 
not allow to resolve the hyperfine components. The file contains 
frequencies, energy levels, and E instein A coefficients from the 
JPL catalog ( Pickett et alj[T998l) and collision rates with ortho- 
and para-H2 for t emperatures in the range of 15 - 300 K from 
Off er et al.1 (11994 . The highest excited state contained in the file 
has an energy of 875 K. The ortho-to-para-H2 ratio in the model 
is temperature dependent according to 



mm 



3.0,9.0xexp 



-170.6 



(1) 



The dust opacities for dust grains with thin ice mantles 
for the wavelength range o f 1 - 1300 fim are taken from 
lOssenkopf & Hennind (11994 Table 1, Col. 5). 

The parameter space explored in the slab models covers 
the full range of physical conditions expected in the proto- 
stellar environment, i.e., temperature range T = 50 - 800 K, 
densities of n = 10^ - 10'^ cm"^, dust column densities of 
Nou = 10"* - 10^^ cm , and OH molecular column densities 
of 10'^ - 10'^ cm"^. Temperatures above 800 K are problem- 
atic because the highest rotational level included in the molec- 
ular data file is 875 K and the collision rates are only avail- 
able up to 300 K. The collisional de-excitation rates are kept 
constant above 300 K, and the excitation rates are calculated 
from the de-excitation rates from the detailed balance relations 
and thus depend on the kinetic temperature through a factor 
exp(-A£'/(^B7'kin))- 
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Fig. 23. Excitation temperature from slab models with Tgas — 200 K and Tdust = 100 K for the OH lines at 79, 84, 119, and 163 fim 
for OH column densities of lO'^, 10l^ and 10'** cm and varying H2 density and dust column density. The dashed line indicates 
the critical density of the upper level of each transition, the dotted line marks where the dust becomes optically thick. 



4.2.2. Slab model results 

A useful quantity to describe the excitation of a molecule is the 
excitation temperature, defined as 

AE 

Te. = T (2) 

where g\ and are the statistical weights of the lower and up- 
per level, respectively, and n\ and their normalized popula- 
tions. The energy of the transition is AE - hv and denotes 
Boltzmann's constant. 

Figure |23] shows the excitation temperature of the OH tran- 
sitions at 79.18 jjm, 84.60 jim, 119.44 yum, and 163.12 yum as a 
function of and A^dust for a model of fixed gas and dust tem- 
peratures (Tcjas = 200 K, Tdust = 100 K) and three values of Aqh- 
The AoH = 10'"* cm"^ model largely represents the optically thin 
case, the Aqh - 10'^ cm"- model the optically thick regime, and 
A^OH = 10'^ cm"^ is an intermediate case. The adopted parame- 
ters are representative for the OH emitting region. 

We notice from Fig. |23]that all transitions behave relatively 
similar with respect to their excitation. Furthermore, this be- 
havior is also comparable for all considered gas temperatures 
(50, 100,200,500, and 800 K). For densities above the critical 
density (~ 10** cm"-'), the excitation temperature approaches 



the kinetic temperature of the gas, as expected, irrespective of 
the dust column density. The populations are thermalized at 
the kinetic temperature. Collisions must therefore be the dom- 
inant excitation mechanism in this regime. For densities below 
the critical density, dust pumping dominates the excitation for 
dust column densities above ~ 10^' cm"^. This is illustrated by 
the fact that Tex approaches Tdust, i-C- the populations are ther- 
malized at the temperature of the dust radiation field. Below 
A^dust ~ 10^' cm"^, both mechanisms contribute to the excitation, 
whereby one order of magnitude increase in «Hi and Adust have 
similar effects on T^^. If the OH column density is enhanced, the 
gas becomes optically thick and line photons are trapped. Then, 
the level populations thermalize at the gas temperature already 
at lower densities because the critical density is reduced with an 
effective critical density of ncr('!"ime) ~ «cr(0)/Tiine. 

In models of equal dust and gas temperature, we find that the 
79 yum cross-ladder transition even experiences supra-thermal 
excitation for low densities (hh, ^ 10^ cm"-'), temperatures 
above 50 K, OH column densities lO'"^ - 10'^ cm"^, and dust 
column densities in the range of 2 x 10^*' - 1 x 10^^ cm"^. 

As for line optical depths, the models show that the 119 yum 
transition reaches the highest optical depth as expected from the 
fact that this transition is connected to the ground state. It is fol- 
lowed by the 84 yum transition, originating from the same rota- 
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Fig. 24. Line flux ratios from slab models with Tg^s = 200 K and Tdust = 100 K for flie OH lines at 79, 84, 1 19, and 163 /im for OH 
column densities of lO''', 10"', and 10'** cnT^ and varying H2 density and dust column density. The shaded area marks the observed 
values. 



tional ladder, just above the 119 fim. The 79 fum line is a cross- 
ladder transition with a much smaller Einstein A coefficient and 
does therefore not become optically thick as easily as the two 
afore mentioned ^113/2 intra ladder transitions. 

Figure |24] presents the line flux ratios for the six possible 
combinations of the observed four transitions. We decided to use 
line ratios instead of fluxes because the ratios are less dependent 
on the source size. The black contours represent the minimum 
and maximum value of the observed ratios and the shaded area 
indicates where the models match the observed values, given in 
Table |5] White areas at high dust column densities indicate re- 
gions where at least one line is in absorption. The observed line 
ratios can be modeled with a wide range of parameters, but the 
observed line ratios still allow us to exclude parts of the param- 
eter space. The very optically thin model (A^oh = 10'^ cm"^) 
fails to reproduce the observed 79 yum/119 fim line flux ratio 
for temperatures below 200 - 500 K and can also hardly match 
the 1 19 /^m/163 fim ratio. Similar difficulties occur for the very 
optically thick case (Aqh = 10"* cm"^), where the observed 
79 ^(m/163 fim ratio can only be found in a very narrow range 
of dust column densities. A direct comparison of the column 
density with observations is however difficult because the model 
does not include any geometrical eff'ects. 

There are basically two main regimes in the intermediate 
case (A^oH - 10'^ cm"^) that are able to reproduce to the ob- 
served values. First, there is the radiatively dominated regime 
at low densities (10^ - 10** cm"^) for dust column densities 
in the range of Adust ~ Ixl0'''-3xl0^' cm"^, slighdy 
depending on temperature (lower Adust for higher temperature 
and vice versa). In the regime where collisions start to become 
important (n > 10^ cm"^) and at low dust column densities 
(A^dust 5 10^'^ cm"-), it is harder to find a set of parameters that is 



able to reproduce all line ratios simultaneously, in particular the 
79 /im/84 fim, 79 /im/163 yum, and 84 /im/163 /^m ratios. For 
densities of 10^ - 10** cm"^, small ranges can however be found, 
especially given the uncertainties of our line ratios. At temper- 
atures below ~ 100 K, a collision dominated solution exists to- 
wards higher densities (« > 10** cm""*) and dust column densities 
above ~ 3 x 10^" cm"^. So given the degeneracies in the models, 
we cannot well distinguish between collisional excitation and 
radiative pumping unless we have additional constraints on the 
density of the emitting gas and the dust column density. 

To constrain which of the two sce narios is more likely, we 
ran additional RADEX slab models dvan der Tak et al.1 120071) 
that allow us to include the actual continuum field of a source 
at a given distance. The line width is chosen to be 10 km s ' and 
the density and temperature at a given distance are taken from 
the co ntinuum model of the source presented in Kristensen et al.1 
(l2012h . The OH column density is a free parameter The contin- 
uum field is only included in the excitation of the OH molecules, 
but not in the line formation, thus assuming a geometry in which 
the OH is not right in front of the line of sight to the contin- 
uum. Models which include the source continuum in both the 
excitation and the line formation yield several lines in absorp- 
tion and cannot reproduce the observations. To study the impor- 
tance of the continuum, models without the continuum field were 
calculated as well. The model fluxes were scaled such that the 
84.60 fim flux is equal to the observed value. 

Figure l25]shows the results for Ser SMMl (n = 8 x 10^ cm"-' 
and T = 90.3 K at a distance of 100 AU) and Nqh = lO'*" cm^^. 
While the modeled fluxes are fairly similar for the transitions 
with an upper level energy less than 300 K, the model including 
the continuum field matches the observed values better than that 
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Fig. 25. Comparison of the rotational diagram for the OH lines 
measured from Ser SMMl (black diamonds) and a RAD EX 
model with the source continuum included in the excitation (red 
plus signs) and a model without the continuum field of the source 
(blue plus signs). The model fluxes were scaled such that they 
match the observed 84.60 fim emission. The model including the 
continuum field reproduces the higher excited lines better than 
the one without. 



lo-"* 

CM 

'g 10"'^ 
10-" 



p = 0.76 



o 



^4 



63 /zm 



p = 0.77 



84 /i,m 



10' 10^ 10" 10' 10^ 10=* 
Continuum (Jy) 

Fig. 26. Correlation of OH 84 yum flux with the continuum fluxes 
measured from the PACS data at 63 ;um and 84 fim. 



without a source continuum for the higher excited lines (98 fim, 
65 fj.m, and 71 /um). 

Another indication that the continuum may be important is 
illustrated by Fig. |26] showing the correlation of the OH 84 fim 
flux with the continuum at 63 fim and 84 fim as measured from 
the PACS data. A continuum field with a temperature similar to 
the observed rotational temperature of OH, i.e. Tiad ~ 100 K, 
would peak at around 30 fim, which cannot be observed with 
Herschel. Therefore, we use the continuum at 63 fim and 84 fim 
instead. 

For the Orion bar, iGoicoechea et al.l (1201 ll) found line ra- 
tios that differ from the values derived for our source sample, 
except for the 79 //m/119 fim and 79 yum/163 ftm ratios. The 
79 /im/84 fim ratio of almost four is higher than in our sam- 
ple, as is the 119 jum/163 fim ratio. On the other hand, the 
84 fim/ 119 fim and 84 jL/m/163 fim ratios are lower The tem- 
perature, density, and the column density of the dust grains are 
lower in the Orion bar, therefore both collisional excitation and 
pumping by the FIR field are likely lower. 



4.3. Discussion 

From the PACS observations we find that the excitation of the 
OH emission in our sample of 17 low- and 6 intermediate-mass 
protostars is similar in all sources, as indicated by the relatively 
constant line ratios among the target sample. This could either 
mean that the OH emission from all sources arises from gas at 
very similar physical conditions or that the observed line ratios 
can be realized by a broad range of gas properties. Comparison 
of the ratios to the modeling results (cf. Sect. 14.2.2b shows that 
the observed values can indeed be reproduced by a wide param- 
eter range, which may explain why similar ratios are observed 
in sources that span a large range in luminosity, envelope mass, 
and evolutionary stage. 

From the full range scans of the two class YSOs Ser SMMl 
and NGC 1333 IRAS 4B we find that the rotational temperature 
of OH in these sources is around 70 K. From the fact that the 
points can be relatively well fitted with a straight line in the rota- 
tional diagram, we conclude that the level population of OH can 
be represented by a Boltzmann distribution at T a; 70 K. The 
constant line ratios found in the sample may be an indication 
that other sources could have similar rotational temperatures. 
However, the models demonstrate that different physical condi- 
tions are able to reproduce the observed line ratios and thus full 
range scan data, which includes more than four OH transitions, 
of additional sources is needed to test how the OH rotational 
temperatures vary from source to source and whether there is any 
difference between evolutionary stages. If the rotational temper- 
ature were equal to the ambient kinetic temperature, the popula- 
tion would be in LTE. However, because of the very high critical 
densities of most OH transitions (around 10** cm""*), this sce- 
nario seems to be rather unlikely, as a significant amount of th e 
emission stems from the outflow (IWampfler et al.ll2010i 1201 ll) . 
Compression factors of a hundred are needed to obtain a post- 
shock density of 10^ cm"^ for a pre-shock density of 10^ cm"-^ 
(Neufeld & Dalgarno 1989a). Alternatively, the OH excitation 
may be in equilibrium with the temperature of the far-infrared 
dust continuum field, that can strongly pump OH. 

The models show that the degeneracy in the parameter space 
is large, allowing different combinations of parameters for the 
realization of the observational results. Therefore, it is only pos- 
sible to distinguish between the collision-dominated and radia- 
tive pumping controlled scenarios described in Sect. 14. 2.2] when 
additional constraints on the density and dust column density 
of gas emitting or absorbing in OH far-infrared lines are avail- 
able. The similar line ratios could be a hint that the excitation 
is not controlled by a combination of both collisions and radia- 
tive pumping but rather by one dominant mechanism, because 
it would be much harder to match the same or similar combina- 
tion in all targets. The RADEX results point towards radiative 
pumping as the main excitation process, at least in Ser SMMl, 
as models including a far-infrared continuum field fit the obser- 
vations better, in particular the high excited lines (65 fim and 
71 fim). Another indication that the radiation field is important 
in the excitation is provided by the strong correlation of the OH 
84 fim luminosity with the bolometric luminosity of the sources 
regardless of their evolutionary stage and type, and a tentative 
trend with the continuum measurements. 

An envelope-only or strongly envelope-dominated scenario 
can be excluded for the low-mass sources from the fact that 
spherically symmetric radiative transfer models yield the 79, 84, 
and 119 fim transitions in absorption, in conflict with the ob- 
served emission lines. Furthermore, the FWHM of the modeled 
163 fim line is at least a factor of two narrower than the ob- 
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served value, pointing also at an outflow origin of the emission. 
Several physical regimes have been identified in the outflow 
from Herschel H2O and CO observations (e.g. Nisini et al. 2010; 
Yildiz et al. 2010; van Kemcen et al. 2010a; Lefloch et al. 2010; 
Kristensen et al., ,2012; .Herczeg et al., ,2012, ; iVasta et al., .2012) : 
the cold entrained outflow material, with typical densities around 
10^ cm ' and temperatures around 100 K, an intermediate warm 
region (T ~ 300 K, n ~ 10*" - 10^ cm^'), and the shock front, 
where temperatures reach T ~ 800 K and densities are found up 
to w ~ 10** cm ^ The OH rotational temperature of T ~ 70 K 
is low compared to the values found for H2O and CO and fa- 
vors the intermediate w arm regime over a shock front origin 
(iGoicoechea et al.ll2012 ). The entrained material however has a 
density that is significantly lower than the critical density, which 
would require the presence of a strong radiation field to excite 
the OH molecules, for instance close to the protostar. 

5. Conclusions 

We analyze HerschelfPACS spectroscopy of the OH lines at 
79, 84, 119 and 163 yum in a set of 23 low- and intermediate-mass 
protostars. In addition, we use slab radiative transfer models to 
study the OH excitation. Our main results are: 

1. At least one OH line is detected in all 23 sources. For 
NGC 1333 IRAS 2 A, only the 1 19 yum doublet is detected in 
absorption against the continuum. Most sources show only 
emission lines. Absorption occurs towards higher envelope 
masses. 

2. The derived rotational temperatures for Ser SMMl and 
NGC 1333 IRAS 4B, where additional OH lines are avail- 
able from a full spectral scan, are similar and around 70 K. 

3. The ratios between the fluxes of different OH emission lines 
are relatively constant among the source sample, i.e. the 
fluxes in the different transitions are well correlated. 

4. There is a strong correlation (Pearson correlation coefficient 
of p = 0.93) between the OH line luminosity and the bolo- 
metric luminosity of the sources. The correlation with enve- 
lope mass is less tight (p = 0.84) and could even be intro- 
duced by an underlying correlation of the bolometric lumi- 
nosity and the envelope mass. We do not find evidence for 
an evolutionary effect, i.e. a strong correlation with the bolo- 
metric temperature or the evolutionary tracer L°|^^/Menv 

5. There are also trends of an increasing OH flux with increas- 
ing [OI] and H2O fluxes. 

6. The OH 84 yum/H20 89 fim ratio is slightly higher in the low- 
mass class I YSOs than in the class sources in our sample. 
The ratio seems to increase with increasing bolometric tem- 
perature and to decrease towards higher envelope masses, in- 
dicating that a change in the ratio could be an evolutionary 
effect. 

7. Spherically symmetric protostellar envelope models do not 
fit the observed OH line widths and produce 79, 84 and 
119 fim lines in absorption rather than emission. 

8. Slab radiative transfer models representative of an outflow 
show that the observed line ratios can be reproduced by a 
range of parameters. The degeneracy in the parameter space 
of the models and the currently available data make it chal- 
lenging to clearly distinguish whether the excitation of OH 
is dominated by collisions or by radiative pumping or even a 
combination of both. Fluxes of transitions with E^p > 300 K 
are clearly better reproduced by RADEX models with ra- 
diative pumping compared to those without a source contin- 
uum. The OH cannot be right along the line of sight to the 



continuum, however, since this would yield several lines in 
absorption. 

9. Both the excitation and spatial extent observed by Herschel, 
combined with independent information on broad line 
widths, provide strong evidence that the OH emission from 
low-mass protostars is dominated by shocks impacting on 
the dense envelope close to the protostar. 
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Table A.l. Observing modes, identities, and pipeline versions for the sources sample. 



Source 


Mode 


Setting 1 


Pipeline 


Setting 2 


Pipeline 


NGC 1333 IRAS 2 A Range 


1342190686 


v6 


1342191149 


v6 


NGC 1333 IRAS 4A Line 


1342216082 


v8 


1342216083 


v8 


NGC 1333 IRAS 4A Range 


1342216084 


v5 


1342216085 


v5 


NGC 1333 IRAS 4B Line 


1342216180 


v8 


1342216179 


v8 


NGC 1333 IRAS 4B Range 


1342216178 


v6 


1342216177 


v6 


L 1527 


Line 


1342225836 


v8 


1342192983 


v8 


CedllOIRS4 


Line 


1342210188 


v8 


1342210187 


v8 


IRAS 15398/B228 


Line 


1342191302 


v8 


1342191301 


v8 


L 483 mm 


Line 


1342217806 


v8 


1342192156 


v8 


Ser SMMl 


Range 


1342207781 


v6 


1342207780 


v6 


Ser SMM3 


Line 


1342207778 


v8 


1342207779 


v8 


L 723 mm 


Line 


1342208919 


v8 


1342208918 


v8 


L1489 


Line 


1342214361 


v8 


1342214360 


v5 


TMR 1 


Line 


1342225835 


v8 


1342225834 


v8 


TMC lA 


Line 


1342225833 


v8 


1342225832 


v8 


TMC 1 


Line 


1342225831 


v8 


1342225830 


v8 


HH 46 


Line 


1342186315 


v8 


1342186316 


v8 


DK Cha (DIGIT) 


Range 


1342188039 


v6 


1342188040 


v6 


RN0 91 


Line 


1342215640 


v8 


1342215639 


v8 


AFGL 490 


Line 


1342202582 


v8 


1342191353 


v8 


NGC 2071 


Line 


1342218760 


v8 


1342218761 


v8 


Vela IRS 17 


Line 


1342209407 


v8 


1342211844 


v8 


Vela IRS 19 


Line 


1342210189 


v8 


1342210190 


v8 


NGC 7129 FIRS 2 


Line 


1342186321 


v8 


1342186322 


v8 


L1641 S3MMS1 


Line 


1342226194 


v8 


1342226195 


v8 



Table A.2. H2O and [OI] fluxes of the intermediate-mass sources in units of 10 W m ^ measured on the full 5x5 array. 



Source H2O transitions [01] transitions 





32,2 " 


-2i,i 


2i,2- 


-lo,i 




■^P2 


^Po- 


■^Pi 




89.99 jum 


179.52 jum 


63.18 jum 


145.53 jum 


AFGL 490 






(abs.) 


121.1 


± 1.4 


11.88 


±0.26 


NGC 2071 


12.80 


± 1.03 


43.70 


±0.65 


552.2 


± 1.2 


42.42 


±0.41 


Vela IRS 17 


2.46 


±0.25 


12.22 


±0.18 


237.6 


±0.79 


56.85 


±0.33 


Vela IRS 19 






3.28 


±0.13 


95.48 


± 1.08 


11.04 


±0.14 


NGC 7129 FIRS 2 


1.84 


±0.38 


4.27 


±0.10 


15.72 


±0.61 


0.90 


±0.13 


LI 641 S3MMS1 


2.71 


±0.24 


1.84 


±0.17 


22.63 


±0.54 


3.03 


±0.12 



Table A.3. Resulting fluxes in units of 10 W m ^ in the spaxel with the strongest continuum (assumed to contain the central source position), corrected by the spillover factor 
(0.70 at 79 /zm, 0.69 at 84 /xm, 0.62 at 119 /xm, and 0.49 at 163 /zm). 



Source 



OH transition [pm] 





79.12 


79.18 


84.42 


84.60 


119.23 


119.44 


163.12 


163.40 


NGC 1333 IRAS 2A (scan) 


- 


- 


- 


- 


abs. 


abs. 


- 


- 


NGC 1333 IRAS 4A 


- 


- 


0.60 ± 0.08 


0.35 ± 0.07 


[0.26 ±0.10] 


(0.37 + 0.10) 


0.26 ± 0.05 


0.56 ± 0.05 


NGC 1333 IRAS 4B 


1.03 ± 0.08 


0.77 ± 0.07 


1.96 ±0.08 


1.18 ±0.07 


1.45 ±0.07 


1.57 ±0.06 


0.33 ± 0.04 


0.35 ± 0.04 


L 1527 


0.33 ± 0.05 


0.27 ± 0.05 


0.50 ± 0.05 


0.49 ± 0.06 


0.60 ± 0.03 


0.71 ± 0.03 


[0.07 ± 0.03] 


(0.12 + 0.03) 


CedllO IRS4 


0.31 ±0.06 


(0.23 ± 0.06) 


0.53 ± 0.06 


0.47 ± 0.06 


(0.16 ± 0.04) 


0.25 ± 0.06 


0.16 ± 0.03 


- 


IRAS 15398/B228 


0.36 ± 0.06 


0.46 ± 0.06 


0.46 ± 0.06 


0.43 ± 0.06 


0.50 ± 0.04 


0.59 ± 0.04 


- 


[0.13 ± 0.05] 


L 483 mm 


1.05 ±0.09 


1.07 ± 0.09 


1.17 ± 0.09 


1.04 ±0.09 


0.78 ± 0.04 


0.95 ± 0.05 


0.25 ± 0.03 


0.22 ± 0.03 


SerSMMl (scan) 


5.71 ±0.69 


blend 


8.57 ± 0.44 


6.25 ± 0.44 


3.70 ± 0.52 


3.29 + 0.51 


2.39 ± 0.34 


1.67 + 0.36 


Ser SMM3 


(0.44 ±0.11) 


(0.36 ±0.11) 


0.63 ± 0.06 


0.42 ± 0.05 










L 723 mm 




[0.17 ± 0.06] 


0.61 ± 0.06 


0.41 ± 0.05 


(0.12 ±0.03) 


(0.13 ± 0.04) 


- 


(0.14 ± 0.03) 


L1489 


0.79 ± 0.07 


0.77 ± 0.06 


1.03 ±0.07 


0.94 ± 0.08 


0.52 ± 0.03 


0.67 ± 0.02 


(0.14 ±0.04) 


(0.15 ±0.04) 


TMR 1 


0.41 ± 0.08 


0.73 ± 0.09 


0.80 ±0.10 


0.88 ±0.10 


(0.13 ±0.04) 


(0.19 ±0.04) 


0.34 ± 0.03 


0.27 ± 0.04 


TMC lA 


0.18 ±0.03 


(0.14 ± 0.03) 


0.72 ± 0.09 


0.71 ± 0.09 


(abs.) 


(abs.) 


(0.10 ±0.03) 


- 


TMC 1 


0.61 ± 0.05 


0.65 ± 0.05 


0.89 ± 0.07 


1.01 ± 0.07 


0.52 ± 0.04 


0.60 ± 0.04 


0.22 ± 0.04 


(0.19 ±0.04) 


HH 46 


0.63 ± 0.09 


0.49 ± 0.09 


0.8 / ± U.08 


0.91 ± 0.09 


V.ZZ + U.U4 


A O /I 1 A A/I 

0.34 ± 0.04 


/A 1 A 1 A AiC\ 

(0.19 + 0.06) 




IRAS 12496/DK Cha (scan) 1.18 ± 0.55 


blend 


2.85 ± 0.34 


2.44 ± 0.30 










RN0 91 


0.61+0.11 


blend 


(0.36 + 0.07) 


(0.40 ± 0.08) 




[0.08 ± 0.03] 






AFGL 490 


(abs.) 


(abs.) 




(abs.) 


abs. 


abs. 


(0.31 ±0.08) 


(0.36 ± 0.08) 


NGC 2071 


10.65 + 0.29 


8.83 ± 0.32 


16.57 ± 0.56 


10.97 ± 0.52 


6.40 + 0.24 


6.96 ± 0.22 


4.09 ±0.17 


3.22 ±0.16 


Vela IRS 17 


1.15 ±0.07 


1.15 ±0.07 


1.18±0.11 


0.61 ±0.10 


abs. 


abs. 






Vela IRS 19 


[0.71 ±0.31] 


(1.26 ±0.26) 


1.75 ±0.40 


1.01 ±0.36 


abs. 


(abs.) 


0.78 ± 0.06 


0.85 ± 0.06 


NGC 7129 FIRS 2 


1.91 ±0.13 


1.56 ±0.14 


3.19 ±0.15 


1.70 ±0.13 






0.55 ± 0.07 


0.67 ± 0.08 


L1641 S3MMS1 


1.78 ±0.18 


1.94 ±0.19 


4.72 ±0.15 


3.38 ±0.15 


(abs.) 


(abs.) 


0.79 ±0.10 


1.17 ±0.10 



CD 

O 



i 

en 
O 

s. 

I 



Table A.4. Resulting fluxes in units of 10 W m ^ in the 3x3 central spaxels. 



Source OH transition [jum] 





79.12 


79.18 


84.42 


84.60 


119.23 


1 19.44 


163.12 


163.40 


NGC 1333 IRAS 2A (scan) 


- 


- 


- 


- 


abs. 


abs. 


- 


- 


NGC 1333 IRAS 4A 


0.73 ±0.13 


(0.64 ±0.13) 


1.72 ±0.15 


(0.48 ±0.14) 




0.94 ±0.13 


0.62 ± 0.04 


0.52 ± 0.04 


JNCjC 1333 IRAS 4ii 


2. i / ± U. 14 


i.o4 ± U.il 


A 1 A 1 
4. 11 ± U.i2 


2.4/ ± U.ll 


2.24 ± U.Uo 


2.78 ± 0.07 


A CI 1 A AT 

0.82+0.07 


A CA 1 A AT 

0.50+0.07 


L 1527 


0.40 ± 0.08 


0.42 ± 0.08 


(0.41 ± 0.12) 




0.40 ± 0.07 


0.52 ± 0.07 




(0.12 ±0.04) 


CedllO 1RS4 


- 


- 


0.77 ± 0.14 


(0.59 ± 0.14) 


- 


- 


- 


- 


IRAS 15398/B228 


0.63 ± 0.14 


0.65 ± 0.14 


1.18 ± 0.13 


0.58 ± 0.12 


0.73 ± 0.09 


0.69 ± 0.08 






L 483 mm 


1.17 ± 0.18 


1.15 ± 0.18 


1.46 ± 0.14 


0.88 ± 0.14 


0.36 ± 0.06 


0.65 ± 0.07 


0.31 ± 0.04 


(0.17 ± 0.04) 


Ser SMMl (scan) 


9.97 ± 0.82 


blend 


12.5 ± 0.43 


9.29 ± 0.43 


8.75 ± 0.71 


8.18 ± 0.71 


2.76 ± 0.30 


1.92 ± 0.33 


Ser SMM3 


1.01 ± 0.18 


1.03 ± 0.17 


1.34 ± 0.18 


1.01 ± 0.17 


(abs.) 


(abs.) 






L, izi mm 






(11.43 ± U.lZ) 


/A A A 1 A 1 A\ 
(U.44 ± U.IU) 










L 1489 


0.82 ± 0.15 


(0.52 ± 0.14) 


0.88 ± 0.12 


0.74 ±0.12 


0.41 ± 0.07 


0.52 ± 0.07 


(0.17 ± 0.04) 


(0.10 ±0.03) 


TMR 1 


(0.72 ±0.16) 


1.17 ± 0.16 


1.35 ± 0.16 


1.12 ± 0.16 


(0.35 ± 0.08) 


0.47 ± 0.08 


0.30 ± 0.05 


0.28 ± 0.05 


TMC lA 




~ 


0.11 ±0.13 


0.66 ± 0.13 


(abs.) 


(abs.) 






TMC 1 


1.00 ± 0.08 


0.89 ± 0.08 


1.40 ± 0.14 


1.11 ± 0.14 


0.45 ± 0.07 


0.67 ± 0.06 


(0.18 ± 0.06) 




HH 46 


0.69 ± 0.17 


0.73 ± 0.18 


0.84 ± 0.17 


1.00 ± 0.18 


0.37 ± 0.06 


0.58 ± 0.06 






IRAS 12496/DK Cha (scan) 




blend 


4.00 ± 0.39 


3.22 ± 0.35 










RN0 91 


(0.52 ±0.17) 


blend 




(0.70 ±0.17) 


(abs.) 








AFGL 490 


abs. 


abs. 






abs. 


abs. 


1.12 ±0.08 


0.91 ±0.08 


NGC 2071 


24.36 ± 0.66 


24.26 ± 0.72 


35.72 ± 0.64 


27.21 ± 0.59 


mix. 


mix. 


10.11 ±0.30 


7.04 ± 0.28 


Vela IRS 17 


1.67 ±0.15 


1.70 ±0.16 


1.81 ±0.22 


1.43 ±0.19 


mix. 


mix. 


0.69 ± 0.08 


0.35 ± 0.07 


Vela IRS 19 


2.50 ± 0.25 


2.46 ± 0.21 


4.71 ± 0.29 


1.88 ±0.26 


mix. 


mix. 


1.13 ±0.06 


0.84 ± 0.06 


NGC 7129 FIRS 2 


2.05 ±0.18 


1.71 ± 0.20 


2.85 ±0.19 


1.57 ±0.16 


0.73 ±0.10 


0.91 ±0.11 


0.71 ± 0.07 


0.51 ±0.07 


L1641 S3MMS1 


1.95 ±0.24 


2.09 ± 0.25 


4.52 ±0.19 


3.13 ±0.19 


mix. 


mix. 


0.79 ± 0.06 


0.93 ± 0.06 



Table A.5. Resulting fluxes in units of 10 W m ^ on the total 5 by 5 detector array. 



Source OH transition [jim] 





79.12 


79.18 


84.42 


84.60 


119.23 


119.44 


163.12 


163.40 


NGC 1333 IRAS 2A (scan) 


- 


- 


- 


- 


abs. 


abs. 


- 


- 


NGC 1333 IRAS 4A 


(0.97 ± 0.23) 


1.22 ± 0.23 


3.11 ± 0.24 


1.17 ± 0.22 


(0.75 ±0.16) 


1.86 ±0.16 


0.75 ± 0.08 


0.63 ± 0.08 


NuC 1333 IRAS 4l3 


2.36 ± 0.30 


1.43 ± 0.25 


5.91 ± 0.35 


'1 on 1 All 

2.87 ± 0.31 


3.17 ± 0.14 


A Art 1 A 1 '1 

4.09 ±0.12 


A A1 1 A 1 A 

0.91 ± 0.10 


/A CO 1 A 1 1 \ 

(0.53 ± 0.11) 


L 1527 










(0.49 ±0.10) 


0.70 ±0.11 






CedllO IRS4 






(0.86 ± 0.27) 




(abs.) 








IRAS 15398/B228 


(1.48 ± 0.34) 


(1.04 ± 0.31) 


1.63 ± 0.29 




1.12 ± 0. 11 


1.25 ± 0.10 






L 483 mm 




(1.22 ± 0.38) 


(0.97 ± 0.28) 




(0.43 ±0.11) 


0.82 ± 0.13 


(0.32 ± 0.07) 


(0.22 ± 0.07) 


Ser SMMl (scan) 


12.20 ± 1.22 


blend 


12.89 ± 0.95 


1.09 ± 0.95 


10.79 ± 0.87 


10.66 ± 0.87 


3.23 ± 0.40 


2.37 ± 0.45 


Ser SMM3 


(1.66 ±0.40) 


2.22 ± 0.38 


2.34 ± 0.38 


1.73 ± 0.36 


abs. 


abs. 






L 723 mm 














(abs.) 




L1489 


(1.06 ± 0.26) 


- 


1.21 ± 0.20 


(0.72 ± 0.20) 


(0.36 ±0.10) 


0.54 ± 0.09 


- 


- 


TMR 1 




- 


(1.02 ± 0.33) 


(1.11 ± 0.33) 


(0.49 ± 0.10) 


(0.45 ± 0.10) 


(0.36 ± 0.08) 


(0.35 ± 0.08) 


TMC lA 


(0.47 ±0.11) 


(0.40 ± 0.12) 


(1.38 ± 0.35) 












TMC 1 


1.32 ± 0.24 




1.44 ± 0.27 


(1.10 ±0.27) 


0.51 ± 0.09 


0.69 ± 0.09 




(0.29 ± 0.09) 


HH 46 


(0.85 ± 0.28) 




(1.02 ± 0.28) 




0.60 ± 0.10 


0.85 ± 0.10 






IRAS 12496/DK Cha (scan) [2.81 ± 1.14] 


blend 


(5.07 ±1.19) 


[2.95 ± 1.06] 










RN0 91 








(1.03 ±0.31) 


(abs.) 








AFGL 490 


mix. 


mix. 






abs. 


abs. 


1.19 ±0.08 


0.90 ± 0.08 


NGC 2071 


30.96 ± 0.73 


30.47 ± 0.79 


42.54 ± 0.81 


33.68 ± 0.75 


mix. 


mix. 


12.62 ± 0.39 


8.33 ± 0.37 


Vela IRS 17 


2.60 ± 0.23 


2.50 ± 0.24 


3.12 ±0.27 


2.92 ± 0.24 


mix. 


mix. 


0.95 ±0.13 


(0.51 ±0.12) 


Vela IRS 19 


3.08 ± 0.44 


2.62 ± 0.37 


4.99 ± 0.49 


2.20 ± 0.44 


mix. 


mix. 


1.27 ±0.07 


1.01 ±0.08 


NGC 7129 FIRS 2 


2.60 ± 0.33 


(1.41 ±0.36) 


1.86 ± 0.30 


(1.30 ±0.26) 


1.05 ±0.11 


1.20 ±0.13 


0.63 ± 0.08 


(0.42 ± 0.08) 


L1641 S3MMS1 


2.49 ± 0.43 


2.93 ± 0.45 


4.28 ±0.31 


2.87 ± 0.30 


mix. 


mix. 


0.87 ± 0.08 


0.93 ± 0.08 
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Appendix B: OH maps 
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Fig. B.l. Maps of the OH emission from NGC 1333 IRAS 2A. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig. B.2. Maps of the OH emission from NGC 1333 IRAS 4A. The red and blue dashed hnes indicate the rest wavelengths of the OH and CO transitions, respectively. Spaxels 
in gray were not included in the flux measurement. 
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Fig. B.3. Maps of the OH emission from NGC 1333 IRAS 4B. The red and blue dashed hnes indicate the rest wavelengths of the OH and CO transitions, respectively. Spaxels 
gray were not included in the flux measurement. 
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Fig. B.4. Maps of the OH emission from L 1527. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig. B.5. Maps of the OH emission from Ced 110 IRS 4. The red and blue dashed hnes indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig. B.6. Maps of the OH emission from IRAS 15398. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig. B.7. Maps of the OH emission from L 483. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig. B.8. Maps of the O H emission from Ser SMMl. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. See also 
iGoicoechea et all ( |2012|) . 
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Fig. B.9. Maps of the OH emission from Ser SMM3. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. Spaxels in gray were 
not included in the flux measurement. 
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Fig. B.IO. Maps of the OH emission from L 723. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig. B.ll. Maps of the OH emission from L 1489. The red and blue dashed hnes indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig.B.12. Maps of the OH emission from TMR 1. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig.B.13. Maps of the OH emission from TMC lA. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig.B.14. Maps of the OH emission from TMC 1. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig.B.15. Maps of the OH emission from HH 46. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig.B.16. Maps of the OH emission from RNO 91. The red and blue dashed hnes indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig. B.18. Maps of the OH emission from NGC 2071. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. 



I 

13 
131 



en 



u 

a 

< 



30 



20 E- 



10 



Vela IRS 17: OH 79//m 



^ '- 



-10 E- 



-20 E- 



-30 



J 



I I II III l| I, " l l - M 'l 



78.9 79.4 



lIj 



l1i 



l1i 



l1i 



30 20 10 -10 -20 -30 
ARa (arcsec) 



u 

cn 
U 

ca 
u 

V 

a 



30 
20 |- 

10 |- 

|- 

-10 |- 

-20 |- 

-30 L 
30 



Vela IRS 17: OH 84/xm 



HI 



T — r 

I I 



1 — r 



m 



T — F- 



T — r 



— 1 — r— 



l1i 



lIi 



lIj 



84.2 84.8 



l1i 



20 10 -10 -20 -30 
ARa (arcsec) 



w 
u 

ee 
u 

0) 

a 

< 



30 



20 '- 



10 E-f 



E- 



-10 E- 



-20 E- 



-30 



Vela IRS 17: OH 119/zm 



Ibi: iJir. „ 



f^" - 111 - llil". 



jir : u ■ I] iri 



J 



J... 



III n 
■'lr'.:".T'" 



.PiilnnnUlj 



,l^l..^J,JI 

" y" uuu m I 



nplLjl 



M 



119.0 119.(1 



v 
n 
u 

CO 

u 

Q 

< 



30 



20 E- 



10 E-f 



30 



20 



-- 

-10 |- 

-20 |- 

-30 L 
30 



Vela IRS 17: OH 163/xm 



162.2 163.e 



10 -10 -20 -30 30 20 10 -10 -20 -30 

ARa (arcsec) ARa (arcsec) 

Fig. B. 19. Maps of the OH emission from Vela IRS 17. The red and blue dashed Unes indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig. B.20. Maps of the OH emission from Vela IRS 19. The red and blue dashed Unes indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig.B.21. Maps of the OH emission from NGC 7129. The red and blue dashed lines indicate the rest wavelengths of the OH and CO transitions, respectively. 
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Fig. B.22. Maps of the OH emission from L 1641 . The red and blue dashed hnes indicate the rest wavelengths of the OH and CO transitions, respectively. 
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